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Abstract. 

Solutions to explicit time-dependent problems in quantum mechanics are rare. 
In fact, all known solutions are coupled to specific properties of the Hamiltonian 
and may be divided into two categories: One class consists of time-dependent 
Hamiltonians which are not higher than quadratic in the position operator, like 
i.e. the driven harmonic oscillator with time-dependent frequency. The second 
class is related to the existence of additional invariants in the Hamiltonian, which 
can be used to map the solution of the time-dependent problem to that of a related 
time-independent one. 

In this article we discuss and develop analytic methods for solving time- 
dependent tunneling problems, which cannot be addressed by using quadratic 
Hamiltonians. Specifically, we give an analytic solution to the problem of 
tunneling from an attractive time-dependent potential which is embedded in a 
long-range repulsive potential. 

Recent progress in atomic physics makes it possible to observe experimentally 
time-dependent phenomena and record the probability distribution over a long 
range of time. Of special interest is the observation of macroscopical quantum- 
tunneling phenomena in Bose-Einstein condensates with time-dependent trapping 
potentials. We apply our model to such a case in the last section. 
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Hamiltonian H(t) 


Reference 


H (t) = & " VoS(x)/t 





H(t) = ^-V S(x)/Vt 2 + Bt + C 


M 


H{t) = £ - \muj 2 x 2 - V e-^5(x) 


this work 



Table 1. Known analytic solutions to time-dependent Hamiltonians H(t) that 
are not quadratic in x,p. 



1. Solutions to time-dependent problems in quantum- mechanics 

In this paper we construct the time-evolution operator (or propagator) K(x, t\x , if) 
for some (time-dependent) potentials in position space. The use of the propagator 
enables us to discuss the time-evolution of an initial state ip(x,t'), which is given at 
time t': 

Surprisingly, it is very hard to find analytic expressions for the propagator (for a 
comprehensive list see pQ). Time-dependent problems in quantum- mechanics are 
divided into two groups: 

(i) The sudden switch from one Hamiltonian to another one at time t', where both 
Hamiltonians are not time-dependent: 

(2) 

The sudden change of the potential will lead to a spatial time-evolution of the 
wave function tp(x,t), since in general a stationary eigenstate of will not be 
a stationary eigenstate of H + . In this sense, H— generates an initial wavepacket, 
which will propagate under the influence of H+ for t > t'. 

(ii) Explicit time-dependence due to a continuously changing Hamiltonian H (t) for 
t > 0. 

In the first case, one can expand the given eigenstate of if_ in terms of the eigenstates 
of H + , or if one models a continuous source of particles (i.e. a particle beam) one can 
use the powerful method of Laplace transforms, either in an approach that directly 
incorporates the initial wave function ip{x,t') 1 or by construction of the energy- 
dependent Green function. Instructive examples are the problem of diffraction in 
time (the Moshinsky shutter) [2] and the evolution of quantum states in external 
fields 011]. 

Unfortunately the method of Laplace transforms does not cope well with time- 
dependent potentials. For this type of problems, the Laplace transform gives rise to 
integro-differential equations. From a conceptual point of view this obstacle is not 
unexpected, since it is no longer possible to work with a conserved energy. Analytic 
expressions for propagators are only known for two remaining subclasses: 

1) quadratic Hamiltonians with time-dependent coefficients (i.e. the driven harmonic 
oscillator with a time-dependent frequency) pQ , 

2) time-dependent potentials that possess a time-dependent scaling parameter (all 
the potentials in Table belong to this class) |HJ |SJ |3| ■ 
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Of special interest for tunneling phenomena are potentials that possess a dip and 
a barrier through which an initially confined quantum state can decay (see Fig. |3J) . 
Commonly used models for tunneling barriers are cubic, Eckart-type, or 6 potentials 

[3 ei ei nni - 

2. Tunneling out of a one-dimensional quantum well 



In this work we will derive a novel analytic solution for the tunneling out of a quantum 
well (modelled by an attractive (5-potential V(x 1 t) = — V(t)S(x)) which is embedded 
into a repelling oscillator potential H n = 2— — ^mui 2 x 2 . We will discuss and compare 
three different cases: 



(i) The propagator for an inverted oscillator potential: 

,2 l 



H = 



P 
2m 



2 2 

—mui x . 
2 



(3) 



(ii) The propagator for a time-independent attractive (5-potential embedded in a 
repelling oscillator: 

- 2 1 

' ' " ' (4) 



P 



H=^-- > 2 i 2 - V S(x). 
2m 2 

(iii) The propagator for an increasingly less attractive (5-potential embedded in a 
repelling oscillator: 

- 2 1 

' ' ' (5) 



H 



"-mu?x 2 - V e-^S{x). 



2m 2 

Since we are interested in the decay rate of the initially confined state, we have 
to define a quantity that will convey this information. One possibility is to calculate 
the ,, survival probability", which is the overlap of an initial state with the state at a 
later time: 



P(t) = 



dx / dx' ip*{x,0)K(x,t\x',0)ip(x',0) 



dxip*(x,0)ip(x,t) 



(6) 



For the initial state, we will always chose the bound eigenstate of the Hamiltonian 

2m 

mV \x\ 



Hs=^~ V S(x) 



VmVo 
Ws{x) = — - — exp 



ft 2 



(7) 



2.1. Propagation under the influence of an inverted oscillator 
The propagator for 



H n = 

is given by setting lo = 
K n {x,t\x',0) 



P 

2m 



1 2 2 

—mui x 
2 



in the propagator for the attractive oscillator ^2 ^ : 
imw (x 2 + x' 2 ) cosh(wi) — 2xx' 



2Trihsmh(ujt) 



exp 



2 sinh(wt) 



(8) 



(9) 
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Figure 1. Survival probability of the initial state %ps{x) = e~\ x \ for different 
potentials: (a) free particle, (b) inverted oscillator and static <5-potential, (c) 
inverted oscillator and exponentially decreasing 5-potential, (d) inverted oscillator 
potential. Parameters: atomic units (m = H = 1), u) = 1, Vq = 1. 



The time evolution of the initial state ips{x) can be expressed in terms of the 
Moshinsky function M(x, k,t) which may be defined in different ways: 

\2~ 



M(x,k,t) 



dx'- 



1 



exp 



— / dn 

2tt 



2nit 
exp [ift 



ikx' 



.(x-x'Y 



2t 



erfc 



tt/4 \ x ~ 

Using equations and (|10|l . wc obtain the wave function at a later time t: 



iexp 



i/cx — 2^ 2 ^ 



(10) 

(11) 
(12) 



m 2 u)Vo 



dx' exp 



imu> (x 2 + x 12 ) cosh(wi) — 2xx' mVo \x'\ 



-M 



mVo 



2irih 3 smh(ujt) J _ 

imVo xmuj 



coth(wi) 



M x 



h 2 



2 sinh(wt) 



1 



n 2 



htSLuh(ujt) 



imVo xmuj 



h 2 



h 



1 



tanh(cji) sinh(wt) J ' mw 
/itanh(wi) \ 



(13) 



tanh(cji) sinh(o;t) / ' muj 

The corresponding survival probability is shown in Fig. ^ Not surprisingly, the 
presence of a repelling potential expedites the decay process compared to the free 
evolution (with u) — > 0) of the wave packet, which is also displayed. 



2.2. Tunneling from a time-independent quantum well 

A more interesting case is the combination of a repelling long-range potential with an 
attractive short-range potential. A general method to derive the propagator for the 
combination of a J-function potential and another potential is based on the solution 
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of the Lippmann-Schwinger equation by using a Laplace transform (for examples and 
a discussion see [3], Section 3.4.). This method is well suited for time-independent 
Hamiltonians of the form 



H = H - V S(x). 
In this case the Lippmann-Schwinger equation becomes 

i 



(14) 



K (x, t\x\ 0) = K o {x,t\x',0) - 



dt" 



dx"K {x, t\x", t") V(x", t") K(x", t"\x', 0) 



K (x,t\x',Q) + ^ / dt"K (x,t\0,t") K(Q,t"\x',0), 



(15) 



where Kg denotes the propagator for the simpler Hamiltonian Hq. In order to solve 
this equation for the propagator K of the full Hamiltonian H , it is of advantage to 
switch to a Laplace transformed picture. The Laplace transform of the propagator is 
conveniently expressed in terms of the energy-dependent Green function G(x,x'] E) 



G(x,x';E) = ilim / dt J&+»,)t/* K (x, t\x', 0). 

lfl ri^O J Q 



(16) 



If we apply the Laplace transform to both sides of equation l|15|) . we can use the 
convolution theorem for a product of two functions (see [TJ, 29.2.8) and obtain the 
Green G function of the full Hamiltonian H in terms of the Green function Go for Ho: 

G (x,0;E)G (0,x';E) 



G(x,x';E) = G (x,x';E)-V - 1 + 
The inverse Laplace-transform yields the propagator: 



K(x,t\x',0) = K (x,t\x',0) 



2tti 



dEe 



-iEt/h 

G (x,0;E)G (0,x'-E) 
l + VoGo(0,0;E) 



(17) 



(18) 



For the case of a <5-function potential on top of an inverted oscillator we have 



H = H + V(x) 



P 



1 



H = ^ -mw 2 i 2 , 

2m 2 



and 



V(x) = -VoS(x). (19) 



The Laplace transform of the propagator for the inverted oscillator gives rise 
to the parabolic cylinder function (we use |15|. 3.547(10) [there is a misprint, the 
correct result is / °° exp{— 2(acothx + px)} sinh 2 ' 5 xdx = ^a l3 T(p — f3)W_ p p + i (4a)] 
and 9.240): 



1 f°° 

G o (x,0;E) = - / dtc lEt ' h K n (x,t\G,Q) 

JO 

2 -3/4+B/(2iriw) / — 



At the origin this expressions simplifies to 



E \ 
2iniu ) ^ «. 



E 



1 1 



2\hw 4' 2 



V2|o 



(20) 



(21) 



where B[s,y] denotes Euler's beta function (see 6.2). In principle, we could now 
carry out the inverse Laplace transform by integrating along the real energy axis. But 
since our main interest is to determine the tunneling rate for this type of potential, 
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another method is available to extract this information. Following the procedure 
suggested by Ludviksson for the combination of a 5-Potential with a linear force field 
|16j . we locate the poles of the denominator in equation (|17|l in the complex energy 
plane. For the inverted oscillator, they are given by the roots of the transcendental 
equation 

1 + V G (0,0;E)=0. (22) 

The determination of the roots (denoted by E^) has to be done numerically. Note, 
that there is another set of poles Er present, which results from the T-function in 
G {x 7 Q-E): 



E^=-^ [ k + -y A: = 0,2,4,.... (23) 

All poles have a negative imaginary part and thus lead to a decaying wave- function for 
t — > oo and fixed locations x, x' . However, none of the poles should be included in an 
evaluation of the inverse Laplace transform, since the original contour runs along the 
real energy axis and cannot be closed in the lower imaginary Emplane: The residues 
at the poles with negative imaginary part correspond to wavefunctions that are not 
bounded at infinity (see also the discussion in {H^)- However, the first pole E^ 
lies closest to the real energy axis and is directly connected to the presence of the 
J-Potential: For u> — > the pole is located at E = — mV§ / (2/i 2 ), which is precisely the 
binding energy of the ^-Potential. The inverted-oscillator potential leads to a shift of 
this pole to 

E^=E r -^Ei, E r , Ei e R. (24) 
An estimate for the lifetime of the decaying state is thus given by 

Tm = (25) 

Fig-fflshows that the exponential decay rate is slowed down compared to the evolution 
in a purely inverted oscillator potential due to the attraction of the ^-potential. 

2.3. Time- dependent 6-potential and an inverted oscillator 

Interestingly, it is possible to obtain an exact analytic expression for the propagator 
of a time-dependent ^-potential on top of an inverted harmonic oscillator. The basic 
idea is to employ a mapping procedure, which transforms the explicit time-dependent 
problem into a related time-independent one. Details of this method are given in 
|H1 Em • The potential must have the following explicit time-dependence: 

V{x ' t) ^ckt) V [ K W) +A{t) ) +^( t ) x + 92(t)x 2 +g (t), (26) 

where the time-dependent coefficients C(t), A(t),go(t), gi(t), g-i (t) are coupled by a set 
of ordinary differential equations: 

C(t)-2g 2 (t)C(t) =K/C 3 (t), (27) 

, 2A{t)C{t) KA(t) 9l (t) 

A{t)+ c(t) + ^W~c(tj- > (28) 

and where K denotes an arbitrary constant. If we choose for the potential the form 
V(x, t) = -^muj 2 x 2 - e- ut V 5(x), (29) 



1 
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we can fulfill all the differential equations by setting 

C(t) = e<"* , A(t)= g (t) = 9l (t) = K = 0, g 2 (t) 



1 2 
-—muj . 

2 



(30) 



Following , we can now construct an equivalent time- independent Hamilton operator 
with a potential V(x) by transforming the position operator and the time according 
to 



= e x, 



x', 



t = 



C(t) C(Q) 
The propagator for the transformed Hamiltonian 



ds 



1 



2u 



P 2 



P 2 



H=^- + V(x) = - V Q S(x) 
Zm Zm 



(31) 



(32) 



function (k — y/2mE/h) 



can be derived with the method of the previous section with the free particle Green 

(33) 



G frcc (x,x';E) = -— e ik \*- x , K hee (x, t\x', 0) 



m 
2-kM, 



exp 



im(x — x') 



2ht 



The inverse Laplace transform is conveniently expressed in terms of the Moshinsky 
function (see i.e. equation (3.43) in 0, replace Vq by — Vq): 



K(x,t\x',0)=K iree (x,t\x',0) 



rriVg 
2nih 2 
mVo 



dk- 



exp 



-ihk-'t 



+ ik(\x\ + \x'\) 



k — imVo/h 2 



K bee (x,t\x',0) + ^-M(\x\ + \x\,+imV a /h 2 ,ht/m). 
nr 



(34) 
(35) 



The mapping is associated with a canonical transformation that gives rise to additional 
factors in the propagator [Hj: 



K(x,t\x',0) 



y/O(t)C(0) 



exp 



2h 



,C(t) 



,(7(0) 





Ljt 




exp 


~~2 


exp 



C(t) ' C(0) 
imuj(x 2 — x' 2 ) 1 



k(x,t\x',0) 



2h 



-2uit 



K e'^x, 



2u 



x'O 



(36) 
(37) 



Inserting equation l|34|) shows that the exponential factor exp [— ^1 is compensated 
by Kf Toe (x, t\x' , 0), and the free propagator is mapped into the one of the inverted 
oscillator (as expected from equation 129|l in the limit Vq — > 0): 

imu)(x 2 — x' 2 ) 1 





Ljt 




exp 


~~2 


exp 



2h 



K hoe (x,t\x',0) = K n (x,t\x',0). (38) 



Therefore we may rewrite equation (|36|l like equation l|18fl as the sum of the propagator 
for the inverted oscillator K n and a contribution due to the time-dependent potential 
Ky. 

K(x,t\x',0) = K n (x,t\x',0)+ K v (x,t\x',0), (39) 

with 



K v [x,t\x',0) 



mV n imu>(x 2 -x' 2 )/(2H)- 

h 2 



~M e 



x + x 



imVb h(l — c 



h 2 ' 



2ujm 



.(40) 



The corresponding survival probability of the initial state {7J) is displayed in Fig. ^ 
As expected, the decreasing ^-potential strength expedites the decay process. The 



Tunneling out of a time- dependent well 



8 




10 15 
x \pm] 



Figure 2. Optical trapping potential for Rubidium atoms in the presence 
of the gravitational force along the rc-axis: V(x) = -Vo exp(-2x 2 /wl) - Fx. 
Parameter: Lower curve Vq = 2400 nK, upper one: Vq = 2286 nK, wq = 27fim, 
F = 103 nK//im. Only the lower curve can support a metastable state. 

asymptotic behaviour for u)t 3> 1 can be extracted by inserting lim^t^oo 5 = 0, and 
lim^t-toot = l/(2w) into equation 



K(x,t\x',0) fn exp[-wt/2] 



3 +imw3; 2 / (2/1) 



- imU3X ' 2 /^K{Q,l/{2w)\x',Q) 



Using equation © we obtain an exponentially decaying survival probability 



P(t) 



dxe" 



dx'e 



(41) 



(42) 



Note that this behaviour is independent of the choice of V(x) in equation (|26l) . 



3. The optical atom laser as a time-dependent problem 

As an application of our method we discuss the time evolution of the quantum-state 
of a Bose-Einstein condensate (BEC) that is confined in a time-dependent trapping 
potential. Recent experiments use optically trapped BECs of Rubidium atoms to 
produce a matter wave of macroscopic dimensions |17) . Optical traps utilize the 
polarizability of atoms to achieve a BEC in the presence of a running optical wave, 
which is generated by a focussed laser beam. To a good approximation, the effective 
optical potential of the laser is given by ^Sl 



V tTa . p (x,y,z,t) = - 



V(t) 



1 + (z/zo)' 



exp 



2{x 2 +y 2 ) 



1 



o i + (z/z y 



where 



TTW 2 



and 



V(t) 



^s-Plaser (^) 



(43) 



(44) 



TTCqCWq 

Typical experimental parameters are reported in |17j : The laser is characterized 
by its minimal beam waist wq (27 /im), its wavelength Ai ascr (10.6 fim) and an 
adjustable power Pi aser (200-0 mW). The polarizability a s of 87 Rb atoms is about 
5.3 x 10~ 39 m 2 C/V. The other constants denote the permittivity of the vacuum eo and 
the speed of light c. 
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Since the laser is operated in the presence of the gravitational field Fe x of the 
earth, the Hamiltonian for atoms trapped in the laser-field becomes 



2(x 2 + y 2 ) 1 



1 + (z/zof 



- Fx, (45) 



H(t) = f ^ ex p 

' 2m l + (z/z ) 2 

where F = rag, (g as 9.81 ms~ 2 , m ~ 87 atomic mass units for 87 Rb). The depth of 
the optical potential V(t) is freely adjustable and can be put as a function of time. If 
one reduces the laser power slowly, one can achieve a beam of atoms that leave the 
trapping potentials once it becomes unstable due to the gravitational force. Figure |21 
shows a cut trough the potential along the x-direction with y = z = for two different 
laser powers. The optical potential forms a quantum well, albeit in three dimensions, 
where the repulsive quadratic potential in (JSJ is replaced by a linear potential. 

In order to model the gradually switching off of the optical trapping potential, 
one should take into account the asymmetry of the potential due to the gravitational 
force. We would also like to use the methods developed in Section |2~51 If we choose 
the potential of the form 

V ch (x) = °£ H , with ch(a, f i,x) = ^(ae^+e-^), (46) 

mch (a, fi, x) * 

we can use the known propagator |19j 
K ch (x, x'; t) = K {iee (x, x', t) + 



2ch(a, /i, x)ch(a, fi, x') 



e i^nt/(2m) „ Ln(x'-x) M (x ~ x', -in, th/m) + e^ x - x ')M(a;' - x, -ip,, th/mfj (47) 

to model a time-dependent potential similar to the one used in eq. I|29() 

V ch (x,t) = -±mw*x 2 + ^WM , with C(i)=e-*. (48) 

A realistic (albeit one-dimensional) picture of the time-evolution of this potential is 
given in Fig. |3J Using the same transformation as before, we can derive the time- 
evolution of an initial state, which we choose in this case to be the bound eigenstate 
in the potential V^x): 



iP ch (x,t = Q)= - ^ ^ (49) 



2ch(a, fj,, x) 

The wave-function at a later time becomes 

/>OC 

i>(x, t) = e~ wt/2 / dx' e ^^ 2 '^' 2 )/(^) Kch f e -«* x 



' x',0) Mx',0). (50) 



2w 

The integration has to be carried out numerically. Fig. [3] shows the resulting time- 
evolution of the probability density for different times. The presence of the potential 
well delays the broadening of the initial wavepacket compared to a force-free evolution 
(also shown in the figure) . This behaviour of the model is not desirable if one wishes to 
model wavepackets of longer spatial extension. The experiments so far tried to create 
an elongated wavepacket. To achieve this goal experimentally requires a careful tuning 
of the time-dependence of the trapping potential |2()| . Our simple one-dimensional 
model shows the coherent control of the initial wavepacket by the changing potential 
and may describe a more focussed atomlaser. 

A one-dimensional treatment as presented here should reproduce the main 
features of the experiment, since the tunneling in the directions perpendicular to 
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Figure 3. Time evolution of an initially confined wavepacket in a time-dependent 
potential. The dotted line shows the potential given in eq. 1481 . the dashed line 
shows the force-free time-evolution of the initial state, and the solid line denotes 
\tp(x,t) 2 \ from eq. JSDJ. Parameter: a = 1/450, (i = 2x 10 6 /m, u) = 2tt x 50 Hz, 
the potential energy scale is nK, and arb. units, for \tf>(x,t)\ 2 . 



the gravitational field are not important, because the lowering of the potential 
is related mainly to the gravitational potential — Fx. For an accurate theory of 
the experimentally realized atomlaser |17| . one should employ a three-dimensional 
trapping potential and a deeper optical potential, which can accommodate several 
bound states. A detailed model of an optical trap by reflect ionless potentials [23 is 
another possibility to obtain potentials which still admit analytic solutions. 

4. Conclusions 

We have conducted a systematic study of tunneling phenomena from an attractive 
potential emebedded into an inverted oscillator barrier. The use of the propagator 
makes it possible to actually obtain the complete time-evolution of the wavefunction 
analytically for an exponentially decaying quantum well. Our results are relevant 
for the study of tunneling processes in time-dependent confining potentials, such as 
the trapping potentials for condensed Bosons (BEC). We succeeded in finding new 
analytic results for one-dimensional tunneling processes through a time-dependent 
barrier. The analytic results provide important reference cases for other methods, 
such as the semiclassical WKB approximation |S]. The mapping procedure of 
Section 12.31 in principle is not limited to one-dimensional problems and opens the 
way to the study the propagation of wavepackets through explicit time-dependent 
barriers. The experimental realization of coherent atomic ensembles in controllable 
confining potentials demonstrates the need for theoretical descriptions of tunneling 
processes through time-dependent barriers. 
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